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Abstract 

Increasing the sensitivity of a gravitational- wave (GW) detector improves our ability to measure 
the characteristics of detected sources. It also increases the number of weak signals that contribute 
to the data. Because GW detectors have nearly all-sky sensitivity, they can be subject to a confusion 
limit: Many sources which cannot be distinguished may be measured simultaneously, defining a 
stochastic noise floor to the sensitivity. For GW detectors operating at present and for their 
planned upgrades, the projected event rate is sufficiently low that we are far from the confusion- 
limited regime. However, some detectors currently under discussion may have large enough reach 
to binary inspiral that they enter the confusion-limited regime. In this paper, we examine the 
binary inspiral confusion limit for terrestrial detectors. We consider a broad range of inspiral rates 
in the literature, several planned advanced gravitational-wave detectors, and the highly advanced 
"Einstein Telescope" design. Though most advanced detectors will not be impacted by this limit, 
the Einstein Telescope with a very low frequency "seismic wall" may be subject to confusion noise. 
At a minimum, careful data analysis will be require to separate signals which will appear confused. 
This result should be borne in mind when designing highly advanced future instruments. 
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I. INTRODUCTION 



Compact binary coalescences, the gravitational- wave (GW) driven inspiral and merger of 
binaries whose members are neutron stars and/or black holes, are among the most promising 
sources of GWs for ground-based detectors. The late inspiral and merger waves from these 
binaries are very strong, and their waveforms are predictable. Their detection would open 
a direct window into strong gravitational dynamics, teaching us much about the evolution 
of massive stellar systems, the nature of gravity in the strong-field regime, and perhaps the 
nature of matter in neutron stars. The well-modeled nature of their waves may allow them 
to be used as standard candles (or standard sirens) [H El El HI E] to probe cosmology. By 
combining GW measurements, which encode distance to the source, with optical or radio 
observations of the host galaxy (for example, if the event is associated with a gamma-ray 
burst), which provides redshift, one can calibrate the relation between luminosity distance 
and redshift (the Hubble diagram) [6]. 

Such science goals have played a large role in driving plans for ever more sensitive GW 
detectors. Improved GW sensitivity impacts the measurement of binary coalescence events in 
two ways. First, any measured event is characterized with greater precision by an instrument 
with high sensitivity thanks to improved signal-to-noise ratio. Many of the analyses of GW 
science from compact binary measurement assume measurement by "advanced" detectors 
(e.g., [HE]), with sensitivity to GWs roughly a factor of ten or so greater than those in 
operation today. Second, the "detection horizon" to which events can be measured grows 
as sensitivity is improved. GW detector noise is normally quoted as a strain sensitivity. 
Doubling an instrument's sensitivity increases its detection horizon by a factor of two, and 
hence the volume to which an instrument is sensitive by a factor of eight. The rate at 
which events are measured can in principle be substantially increased by relatively modest 
improvements in instrumental sensitivity. 

Such improvement in sensitivity is on the one hand mandated by the rarity of strong 
gravitational-wave events. Extrapolation from the observed binary pulsar systems in our 
galaxy and population synthesis calculations each suggest that binary mergers occur roughly 
once every hundred thousand years per Milky Way equivalent galaxy; a more precise dis- 
cussion of the merger rate is given in Sec. [IT] below. Combining this rate with the density 
of galaxies in our universe drives us to the need for a detection horizon of a few hundred 
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Mpc in order to measure multiple NS-NS events per year. On the other hand, there can be 
too much of a good thing: If the detection horizon is sufficiently far out, sources may be 
detected so often that they overlap, making a transition to a confusion-limited background. 
The transition from discrete sources to confused background has been discussed extensively, 
especially in the context of sources for the space-based detector LISA [9} [TO] . 

In this paper, we investigate the threshold redshift at which merging binary sources for 
ground-based detectors begin to overlap and result in a confusion background. Beyond 
this redshift, it will become difficult to resolve sources individually. As has been found 
in the LISA Mock Data Challenges, separation may be possible by globally fitting for all 
sources simultaneously (see, e.g., Ref. [TT]): similar analyses have been done for the proposed 
highly-advanced Big Bang Observatory [12]. Careful analysis will be needed to see how well 
such techniques can be applied to ground-based detector data. At any rate, it is likely 
that improving sensitivity to probe beyond this point will not be worthwhile. We first 
discuss estimates of the cosmic coalescence rates of neutron star-neutron star (NS-NS) and 



neutron star-black hole (NS-BH) systems in Sec.jlTJ In Sec. Ill we then discuss the different 
detection regimes which pertain to binary coalescence measurement. Key to this discussion 
is understanding how long different sources are in band, which depends most strongly on the 
low-frequency sensitivity of the instrument being used to measure the waves. We consider the 
initial and planned advanced LIGO and Virgo detectors, and three possible configurations of 



the Einstein Telescope concept. Our conclusions are given in Sec. IV Chief among them is 
that confusion issues are likely to become important for ground-based GW detectors which 
are sensitive to these events to a redshift z ~ 1, especially if such instruments have a "seismic 
wall" at or below roughly 5 Hz. This may have important ramifications for designing highly 
sensitive future detectors such as the Einstein Telescope. At a minimum, it will be necessary 
to carefully design the data analysis to disentangle these potentially confused sources. 



II. COALESCENCE RATES 



To motivate the choices for our calculations, we begin by reviewing our current under- 
standing of the rate of NS-NS and NS-BH coalescences in the universe. The final merger 
of a compact binary occurs after two massive stars in a binary system have collapsed to 
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form neutron stars or black holes 1 and have inspiralled through the emission of gravitational 
waves. We assume that the coalescence rate tracks the star formation rate, albeit with some 
delay from formation of the binary to final merger. Accordingly, we put 

m = m x . (i) 

In Eq. Q, p°(z) represents the rate at which binary systems are observed to merge at 
redshift z, and p°(0) is this rate in our local universe. This normalization reproduces the 
local rate for z — 0. We review estimates of p°(0) in more detail later in this section. The 
quantity p* )C relates the past star formation rate (SFR) to the rate of binary merger. It is 
defined as 

In Eq. @), p\ is the SFR, measured in M Mpc 3 yr 1 . The redshift z describes when 
our compact binary merges, and Zf is the redshift at which its progenitor binary formed. 
These redshifts are connected by the time delay t^, which is the sum of the time from initial 
binary formation to evolution into compact binary, plus the merging time r m by emission of 
gravitational waves. It is also the difference in lookback times between Zf and z: 

1 f Zf dz' 

t d = tUzf)-tUz) = W J z {1 + z , )mz ,y (3) 

where 

E(n,z) = ^n A + n m (i + z) 3 ■ (4) 

We use the 737 cosmology [13], with Q m = 0.3, = 0.7 and Hubble parameter H = 70 
km s" 1 Mpc -1 (or equivalently h = 0.7). This corresponds to the concordance model 
derived from observations of distant type la supernovae [13] and the power spectrum of the 
cosmic microwave background fluctuations [13]. Finally, in Eq. ^ P{td) is the probability 
distribution for the delay td- The factor 1/(1 + zj) accounts for time dilatation due to the 
cosmic expansion. 

We find the merger rate per unit redshift by multiplying p°{z) with the element of co- 
moving volume: 

f(^MfM, (5) 



1 We neglect the possible production of compact binaries through interactions in dense star systems. 
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FIG. 1: Cosmic star formation rates (in M Mpc 3 yr 1 ) used in this paper: Hopkins and Beacom 
2006 [19J (continuous line), Fardal et al. 2007 (20] (dashed line), Wilkins et al. 2008 (dot-dashed 
line), |21j . and the fossil model of Nagamine et al. 2006 (dot line). As discussed in the text, these 
rates are largely the same up to z ~ 1, but show important differences at higher redshift. 

where 

dV , . A c r(z) 2 

Note that Eq. Q assumes spatial flatness, Vt m + VL\ = I. 

For p*, we use the recent SFR of [19J, which is derived from new measurements of the 
galaxy luminosity function in the UV (SDSS, GALEX, COMB017) and FIR wavelengths 
(Spitzer Space Telescope), and is normalized by the Super Kamiokande limit on the electron 
antineutrino flux from past core-collapse supernovas. This model is expected to be quite 
accurate up to z ~ 2, with very tight constraints at redshifts z < 1 (to within 30 — 50%). To 
account for uncertainty in reconstructing the SFR, we also considered the model described 
in Fardal et al. (Ref. [2D])- That work uses a different set of measurements and a different 
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dust extinction correction. The SFR found in [20] is the same as that of [19] up to z ~ 1, 
but decreases slightly at higher redshifts. We also consider the model described by Wilkins 
et al. in Ref. [2T] , which is derived from measurements of the stellar mass density. The SFR 
is equivalent to that in [T9J, [20] for z < 0.7, but again is lower at higher redshifts. Finally, 
we consider the SFR of Ref. [22], which is derived from the fossil record of star formation 
in nearby galaxies. It is probably underestimated at small redshifts, and is constant at 
high redshifts due to the contribution of elliptical galaxies. Note that at present there is a 
discrepancy between the "instantaneous" SFR, measured from the emission of young stars 
in star forming regions, and the SFR as determined from extragalactic background light. 
This could have an important impact on the contribution to the confusion background for 
sources from z > 2. 

Population synthesis [21 ESI ESI EH EH E3 [30] suggests that the delay time t d is well 
described by a probability distribution of the form 

Pd(td) oc — with t d > t (8) 

for some minimal delay time r . This broad model accounts for the wide range of merger 
times observed in binary pulsars; it is also consistent with short gamma ray burst observa- 
tions in both late and early type galaxies [21]. Following [3"0~] 152] , who identify a channel 
which produces tight NS-NS binaries with merger time in the range r m ~ 0.001 — 0.1 Myr, 
we assume a minimal delay time for NS-NS of tq ~ 20 Myr. This corresponds roughly to 
the time it takes for massive binaries to evolve into two neutron stars. For NS-BHs, we take 
a minimal delay of r ~ 100 Myr [33], corresponding to the wider orbits and longer merger 
times predicted by classical channels. 

The local merger rate per unit volume, p°(0), is usually extrapolated by multiplying 
the rate in the Milky Way (r mw ) with the density of Milky- Way equivalent galaxies. That 
density, in turn, is found by measurements of the blue stellar luminosity to be roughly 
n mw ~ (1 — 2) x 10~ 2 Mpc~ 3 [311 [351 [36] . Current estimates of the NS-NS galactic coalescence 
rate which extrapolate from observed galactic NS-NS find a rate in the range 17 — 292 Myr -1 
(95% confidence interval), with a most likely rate of 83 Myr -1 |37j. Population synthesis 
predicts a NS-NS merger rate in the range 1 — 300 Myr -1 (most likely near 10 — 30 Myr -1 ; 
see Table [T] for a summary). Rate estimates for NS-BH systems are one or two orders of 
magnitudes smaller, ranging between 0.1-30 Myr -1 (Tab. [J. 
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TABLE I: Current estimates of the Galactic merger rates of NS-NSs and NS-BHs in units of 
Myr" 1 , adapted from Table 4 of [38]. The rates from Ref. [37] are derived from statistical studies; 
the numbers in parenthesis give the 95% confidence limits around the maximum likelihood value. 
All other values are estimated using population synthesis. The high rate obtained in Ref. [SU] is 
due to the assumption that neutron stars or black holes are born with no kick velocity, leading to 
an overestimate of the number of systems that survive both supernovae. The low rate obtained by 
Ref. [13] is due to their treatment of common envelope binding physics. 



Statistics 

Kalogera et al. (2004) [37] 


NS-NS 
83 (17 - 292) 




Population synthesis 


NS-NS 


NS-BH 


Tutunov and Yungelson (1993) |39] 


300 


20 


Lipunov et al. (1997) 00] 


30 


2 


Portegies Zwart and Yungelson (1998) |41| 


20 


2 


Nelemans et al. (2001) [12] 


20 


4 


Voss and Tauris (2003) g3] 


2 


0.6 


de Freitas Pacheco et al. (2006) [27] 


17 




Belczinsky et al. (2007) [S] 


10-15 


0.1 


O'Shaughnessy et al. (2008) [26] 


30 


3 



It is worth noting the stellar mass fraction which produces massive binaries is expected 
to be larger in early-type galaxies than in spiral galaxies thanks to their natter initial mass 
function. Despite their absence of recent star formation, early-type galaxies may make a 
larger contribution to the merger rate than spiral galaxies due to binaries with very long 
coalescence times which were born in a galaxy's first 1 — 2 Gyr. Assuming that ellipticals 
represent a fraction of about 40% of the galaxies, [27] find a small correction factor ~ 2. 

Taking these uncertainties into account, we consider local rates in the range p° c (z = 0) = 
(0.01 - 10) Myr- 1 Mpc" 3 for NS-NS and p° c (z = 0) = (0.001 - 1) Myr" 1 Mpc" 3 for NS-BH. 
Our two reference models for NS-NS are 

• p°(z = 0) ~ 1 Myr -1 Mpc -3 , corresponding to the most probable galactic rate of 
^mw = 83Myr _1 , and 
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FIG. 2: Cosmic coalescence rate, normalized to the local value p°{z = 0), for models with a 
distribution of the delay of the form P(td) oc 1/td- The continuous line is the distribution for a 
minimal delay tq = 20 Myr, which is assumed to be representative for NS-NS binaries. The dashed 
is for a delay of 100 Myr, taken to be representative of NS-BHs. In both cases, we have assumed 
the cosmic star formation rate of Hopkins and Beacom [TH] . 

• = 0) ~ 0.4 Myr -1 Mpc -3 , corresponding to the prediction of the latest population 
synthesis of [26J (r mw = 30 Myr -1 or to r mw = 15 Myr -1 [27J EI]> with a correction 
factor of 2 due to ellipticals). 

Our reference model for NS-BH is p°(z = 0) ~ 0.04 Myr -1 Mpc -3 , corresponding to the 
most recent predictions of [51 J. For all three models, we assume n mw = 1.2 x 10 -2 Mpc -3 

m- 

III. DETECTION REGIMES 

Turn now to the measurement of signals from NS-NS and NS-BH binaries. The contri- 
bution of these binaries to the instrumental data falls into three statistically very different 
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regimes, depending mostly on the typical interval between events (see [52j [53] and references 
within): 



1. Shot noise: This case describes when the number of sources is small enough that the 
interval between events is long compared to an individual event's duration. Measured 
waves are separated by long stretches of silence and can be resolved individually. This 
case pertains to instruments that are only sensitive to events at low redshift. 

2. Popcorn noise: As the reach of instruments increases, the time interval between events 
may come closer to the duration of single bursts. Events may sometimes overlap, 
making it difficult to distinguish between them. 

3. Continuous: For instruments with very large reach and excellent low frequency sensi- 
tivity, the interval between events can be small compared to the duration of an event. 
The signals overlap to create a confusion noise of unresolved sources. 

The average number of measured events at a given moment is given by the duty cycle A. 
It is defined as the ratio, in the observer frame, of the typical duration of a single event r, 
to the average time interval between successive events 

r(z')^(z')dz' . (9) 

Here, dR°/dz' oc p° is the coalescence rate per unit redshift, given by Eq. ([5j; r(z') is the 
observed duration of a GW signal generated at redshift z' . This duration is given at leading 
order by 

= 2567r 5 8/3 G 5/3 [( 1 + ^m C ]- 5/3 fL m , (10) 

where m c = (mim2) 3//5 /( m i + m 2 ) 1//5 is the binary's chirp mass, and where /l is the lower 
frequency bound of the detector, assumed to be much smaller than the frequency at the 
time of the merger. 

The signal duration depends very strongly on both chirp mass and lower frequency bound; 
in particular, signals measured by instruments with very low fi may have very long dura- 
tions. We will take typical NS masses to be 1.4 M , and typical BH mass to be 9.5 M , 
yielding chirp masses m c = 1.2 M for NS-NS and m c = 2.9 M for NS-BH. These values 
agree well with the most probable values derived by jH] with the Star Track population 
synthesis code, for their reference model. The lower frequency bound is determined by 



the properties of the detector used for the measurement. Sensitivity curves describing the 
instruments we include in our analysis are shown in Fig. [3j We consider: 



• Initial LIGO: The present LIGO interferometers [35] have a low-frequency seismic 
"wall" at roughly 40 Hz. NS-NS binaries are in band for roughly 25 seconds; NS-BH 
binaries are in band for roughly 5.8 seconds. The isotropic detection horizon (the angle- 
averaged distance to which a binary can be measured, defined carefully below) is at a 
luminosity distance of 15 Mpc (z = 0.0035) for NS-NS, 30 Mpc for NS-BH (z = 0.007) 
for the LIGO Hanford 4 kilometer detector. When both 4 kilometer detectors (at 
Hanford, Washington and Livingston, Louisiana) and the single 2 kilometer detector 
(at Hanford) are combined (assuming independence of all detectors and stationary 
noise), these numbers increase by a factor of about a/1 2 + l 2 + (1/2) 2 = 1.5 to 22 
Mpc (z = 0.005) for NS-NS, 44 Mpc for NS-BH (z = 0.01). 

• Virgo: The Virgo design sensitivity [46] is comparable to the LIGO sensitivity, but 
with advanced seismic isolation provided by the so-called "superattenuator" (see, e.g., 
Ref. [50] for discussion). A low-frequency wall at 10 Hz means that NS-NS signals will 
be in band for 16.7 minutes; NS-BH signals last for 3.9 minutes. The Virgo sensitivity 
is slightly poorer than the LIGO sensitivity between 100 — 300 Hz, so its isotropic 
detection horizon is 11 Mpc (z = 0.0025) for NS-NS, 23 Mpc for NS-BH (z = 0.0055). 

• Advanced LIGO /VIRGO: The advanced LIGO design [37J |38] should move its wall 
down to a frequency of roughly 10 Hz. NS-NS signals will be in band for about 16.7 
minutes and have a single detector detection horizon of 200 Mpc (z = 0.045); NS- 
BH signals will be in band for 3.9 minutes, and have a detection horizon of 420 Mpc 
(z = 0.09). The current Advanced LIGO design plan will have three 4 kilometer 
instruments; combining their data brings the isotropic horizon to 355 Mpc (z = 0.08) 
for NS-NS and 765 Mpc (z = 0.16) for NS-BH. The plans for advanced Virgo sensitivity 
[4*9"] are comparable to those for LIGO. The expected isotropic detection horizon for 
advanced Virgo is at 150 Mpc (z = 0.035) for NS-NS, 310 Mpc for NS-BH (z = 0.07). 

• The Einstein Telescope: Several possible Einstein Telescope designs are presently un- 
der discussion. The plan is to move the seismic wall down to a frequency below 5 
Hz. A planned low-frequency wall at 5 Hz would mean that NS-NS signals will be in 
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FIG. 3: Strain sensitivities, in Hz" 1 / 2 , for the instruments that we consider as functions of fre- 
quency: LIGO 4 km as specified in the Science Requirement Document |45| : design initial Virgo 
|4"6] : planned Advanced LIGO and Virgo (33 EEJ 09] (optimized for NS-NS detection); and a 
possible Einstein Telescope design (L-shaped, 10 km) [ET science team, private communication]. 

band for 1.8 hours; NS-BH signals will be for 24.6 minutes. Another possibility will 
be to lower it to 3 Hz (so that NS-NS signals are in band for 6.9 hours, NS-BH for 
1.6 hours), or even to 1 Hz (NS-NS lasting for 5.4 days, NS-BH for 1.2 days). The 
Einstein Telescope detection horizon for NS-NS signals is at z ~ 1; for NS-BH, it is at 
z ~ 2. 

The isotropic detection horizon for these instruments is the distance at which the angle- 
averaged signal-to-noise ratio p, defined by 

\kf)\ 2 




KU) 



:df , (11) 



equals the threshold value p t h = 8. In Eq. (11), h(f) is the Fourier transform of the inspiral 



waveform, h n (f) is the strain noise (in units of Hz 1//2 , as shown in Fig. [3J), and () denotes 
an appropriate angle averaging; see, for example, Ref. [6], Sec. II for discussion. We are 
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NS-NS (f =10 Hz) 



10 



10 




FIG. 4: Duty cycle as a function of redshift for NS-NSs (top) and NS-BHs (bottom), for a lower 
frequency bound of 10 Hz and for our reference SFR [19] . The continuous black and grey lines cor- 
respond to the most optimistic and pessimistic estimates of the local coalescence rate, respectively; 
the dashed and dot-dash lines correspond to our reference models (cf. Sec.|ll]). The horizontal line 
at A(z) = 1 indicates the transition between resolved sources and a "popcorn" background; the 
line at A(z) = 10 is our more conservative transition to a confused background. 



careful to specify the "isotropic" detection horizon in our discussion (in reference to this 
angle averaging) to contrast with the term "detection horizon." This name is typically used 
to denote the distance to which an optimally oriented and position source could be detected. 
Whereas fortunate orientation and position may make a source visible beyond the isotropic 
detection horizon, no source is visible beyond the detection horizon. The threshold p t h — 8 
is chosen to keep the false alarm rate (for Gaussian noise statistics) acceptably low. 

The redshift at which sources start to overlap to produce a confusion background is 
the transition between the shot noise and the popcorn noise regimes. This is defined by the 
condition 

A(*,) = 1 . (12) 

In this regime, where the number of sources that overlap is small, it may still be possible to 
separate the sources individually with adequate data analysis techniques. We therefore con- 
sider in addition a more conservative threshold z^, corresponding to the transition between 
the popcorn regime and the Gaussian stochastic background. It is defined by A(z**) = 10. 
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10 



NS-NS (f L =5 Hz) 



10 



10 




FIG. 5: Same as Fig. [4] for a lower frequency bound of 5 Hz. 





FIG. 6: Same as Fig. [4] for a lower frequency bound of 1 Hz. 



Figure [4] shows A(z) for different estimates of the local coalescence rate of NS-NSs and 
NS-BHs and for an instrumental seismic wall of 10 Hz; Figs. [5] and [6] show the same thing 
for seismic walls at 5 and 1 Hz, respectively. The threshold redshifts z* and z„ for NS-NS 
are shown in Table [IT] Table |TTT| shows these thresholds for NS-BH. For ease of comparison, 
we repeat in these tables the detection horizons for the various instruments discussed above. 



IV. DISCUSSION AND CONCLUSIONS 



The main concern of the first and second generations of interferometric gravitational-wave 
detectors will be improving sensitivity to ensure unambiguous first detection of gravitational 
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TABLE II: Top entries: threshold between resolved and unresolved NS-NS binaries for different 
estimates of the source rate p° and detector lower frequency bound /l- No value means that the 
duty cycle is always < 1 or < 10. Bottom entries: lower frequency sensitivity bound and detection 
horizon for current and planned gravitational-wave detectors. For LIGO, the horizon is given for 
the SRD 4 km detector and the value in parenthesis corresponds to the combination of the three 
detectors. 



h 


P°c 








10 


0.01 










0.4 




0.8-0.9 






1 




0.5-0.6 


> 2 




10 




0.2 


0.5-0.6 


5 


0.01 




- 


- 




0.4 




0.4 


1-1.2 




1 




0.25 


0.6-0.7 




10 




0.1 


0.25 


1 


0.01 




0.3 


0.8 




0.4 




0.08 


0.2 




1 




0.06 


0.13 




10 




0.03 


0.06 


LIGO: 




h = 40 Hz 




z DH = 0.0035(0.005) 


Virgo: 




h = 10 Hz 




z DU = 0.0025 


Advanced LIGO: 




Sl = 10 Hz 




zdh = 0.045(0.08) 


Advanced Virgo: 




h = 5 Hz 




zdh = 0.035 


Einstein Telescope: 




/i = {5,3,l} Hz 




Z D H ^ 1 



waves, followed by measurement with sufficient precision to inaugurate gravitational-wave 
astronomy as an observational science. As the technology improves and these instruments' 
sensitivity increases, the possibility emerges that a confusion background may become an 
important effect, ultimately limiting the capabilities of these instruments. 

Referring to Tables [TT| and |III[ we see that there is little danger of a confusion background 
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TABLE III: Same as Table. [Hi but for NS-BH binaries. 



f L Pc 5 ^* ^* * 

10 0.001 



0.04 

1 1.1-1.4 



5 0.001 

0.04 

1 0.5 > 1.6 

1 0.001 > 2.3 

0.04 0.3 0.8-0.9 

1 0.1 0.2 

LIGO: f L = 40 Hz z DH = 0.007(0.01) 

Virgo: f L = 10 Hz z DH = 0.0055 

Advanced LIGO: f L = 10 Hz z DH = 0.09(0.16) 

Advanced Virgo: fi = 5 Hz zdh = 0.07 

Einstein Telescope: fi = {5, 3, 1} Hz zdh — 2 

impacting these instruments even into the advanced LIGO era. The confusion background 
only has an impact when a detector's "reach" (summarized by the detection horizon £dh) 
becomes large and its low frequency "wall" (the cutoff frequency fi) becomes low. A distant 
detection horizon can greatly increase the number of sources that an instrument can reach; a 
low frequency "wall" can greatly increase the amount of time that a source spends in band. 

The combination of these two effects is summarized by the duty cycle, A(z), which is 
the ratio of the typical duration of a measured event to the typical interval between events. 
We find that, for advanced detectors with excellent low frequency sensitivity (such as the 
planned advanced Virgo design), the redshift z* at which the duty cycle for NS-NS is unity, 
A(z*) = 1, may occur close to the detection horizon for the reference coalescence rate. The 



redshift z* defines a "popcorn" background in the language of Sec. Ill The conservative 
redshift z** for which A(z**) = 10 (defining a Gaussian stochastic confusion background) 
may even be close to the detection horizon for optimistic event rates. 
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We find that, for NS-NS, both z* and Zip* oX6 well within the horizon of the planned 
Einstein Telescope, if its low frequency sensitivity is at fz — 1 Hz (ET1). If f L = 5 Hz 
(ET5), we expect the popcorn background to occur before the detection horizon, and more 
likely around ~ 0.25 — 0.4, unless our most pessimistic coalescence rates are accurate 
(p° < 0.015 Myr^ 1 Mpc~ 3 ). The transition to a Gaussian stochastic most likely occurs at 
2** ~ 0.6 — 1.2, but can fall beyond the detection horizon if p° < 0.15 Myr -1 Mpc~ 3 . 

Our conclusions for NS-BH binaries are similar to those for NS-NS for ET1, z* and z** are 
both more likely to occur well below the horizon. For ET5 however there is not likely to be 
enough sources to create a Gaussian stochastic background (even a popcorn background), 
except for the most optimistic coalescence rates (p° > 0.6 Myr -1 Mpc -3 ). As a consequence, 
NS-BHs may always be resolved provided that we can separate them with adequate data 
analysis strategies in the popcorn regime. 

This result motivates very careful analysis of how data from an Einstein- Telescope-type 
instrument would be analyzed. Experience from the Mock LISA Data Challenges [ill] and 
ideas developed for the Big Bang Observatory [12] prove that disentangling multiple signals 
in a gravitational- wave detector's datastream is certainly possible. However, beyond the 
proof of the concept, it is not clear how many of the lessons from these examples carry over 
to the case of ground-based detectors. In particular, even for an Einstein- Telescope-type 
detector, the bulk of the potentially confused signals would be of low signal-to-noise ratio 
since most events will come from the volume of the universe near the detection horizon. It 
would be a valuable exercise to examine how well techniques that simultaneously fit multiple 
signals do in the ET regime. 
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